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Abstract 

In order to provide a microscopic description of levels and E2 transitions in rapidly rotating 
nuclei with internal excitation energy up to a few MeV, use is made of a shell model which 
combines the cranked Nilsson mean-field and the residual surface delta two-body force. The 
damping of collective rotational motion is investigated in the case of a typical rare-earth nucleus, 
namely ^®®Yb . It is found that rotational damping sets in at around 0.8 MeV above the yrast 
line, and the levels which form rotational band structures are thus limited. We predict at a 
given rotational frequency existence of about 30 rotational bands of various lengths, in overall 
agreement with the experimental findings. The onset of the rotational damping proceeds quite 
gradually as a function of the internal excitation energy. The transition region extends up to 
around 2 MeV above yrast and it is characterized by the presence of scars of discrete rotational 
bands which extend over few spin values and stand out among the damped transitions, and 
by a two-component profile in the Ej — E~, correlation. The important role played by the 
high-multipole components of the two-body residual interaction is emphasized. 

PACS-. 21.10.Re, 21.10.-n, 23.20.Lv 
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1 Introduction 

A heavy ion fusion reaction produces a hot and rapidly rotating compound nucleus. After emission 
of neutrons has cooled down the system, the compound nucleus still keeps high angular momenta 
(tens of h) and moderate heat energy, that is a moderate excitation energy (a few MeV) once the 
energy of rotational motion has been subtracted. It then emits a multitude of y-rays, decreasing its 
spin and heat energy gradually, reaching finally the ground state. 

For nuclei with mass number A ~ 170 whose ground state has a stable elongated shape (/3 ~ 0.3), 
up to around 20 different rotational bands have been observed in discrete gamma ray spectroscopy 
studies. The observed rotational bands lie in the lowest excitation energy region near the yrast 
line, representing the “cold” part of the whole of nuclear levels. In the excited region, say at 1 to 

2 MeV above the yrast line, the level density at a given spin is as large as 10^ to 10^ MeV“^ in 
normally deformed rare-earth nuclei. The equivalent level spacing of 10 to 1 keV is rather small 
in comparison to the typical size lOkeV) of matrix elements of the residual two-body nuclear 
force. The residual interaction then becomes effective to cause mixing of many-particle many-hole 
(or many-quasiparticle) configurations in the rotating deformed mean-field potential. Since different 
configurations respond differently to the Coriolis force due to different single-particle alignments, the 
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configuration mixing results in a dispersion of the rotational frequency within each energy eigenstate, 
implying a damping of the collective rotational motion §1- The gamma-rays which are emitted 
from the “warm” region (from about 1 MeV to about a few MeV above yrast, corresponding to 
temperature of a few hundred keV ) cannot be distinguished as discrete peaks and form a quasi¬ 
continuum in the spectra. Rotational damping has been studied experimentally through the analysis 
of the quasi-continuum spectra . Recently, new experimental techniques have been devised 

to deal with the multi-dimensional quasi-continuum gamma-spectra obtained by double or triple 
coincidence gamma-ray experiments, making it possible to study in detail various aspects of the 
collective rotational motion in the warm region I- a 11 B pd] |. In particular, the fluctuation 
analysis method §0 has revealed that the number of rotational bands existing in a rare-earth 
nucleus is only around 30 at a given rotational frequency, thus confirming the occurrence of the 
rotational damping. 

Early theoretical studies of the rotational damping dealt with the E2 strength function associated 


that the configuration mixing 
1,1 li, 111- In this paper we 


with the damped rotational motion, which is obtained by assuming 
is described by the general statistical theory of random matrices [|] 
present a microscopic shell model able to describe individual nuclear levels and E2 transitions in the 
warm region, in order to study in more detail the transition from the regime of discrete rotational 
bands into the regime of damped rotational motion. Our description is based on a mean field 
calculated by the cranking model, which has been quite successful in the description of the low-lying 
portion of the nuclear spectrum (cf. e.g. 0 lii, n, 0 0). We then add a two-body effective 
interaction, mixing the mean-field configurations. The first attempt along this line was made in 
Ref. 10, where a rather schematic residual interaction was employed in the calculations. We instead 
adopt somewhat more ’’realistic” effective nuclear force, i.e., the surface delta force |0|. We found 
in previous studies pl|, |0] that the high-multipole components, which are present in this force, are 
essential to produce the rotational damping. In the present formulation, the emphasis is put not 
only on the residual two-body force but also on the construction of the diabatic single-particle basis, 
which enables us to describe long sequences of the rotational E2 transitions. To be noted is that 
both the particle-rotor model |0| and the interacting boson model |24 have been used to attack 
problems similar to the one we are dealing with in the present paper. Both of these models are more 
complete than the present one with respect to the angular momentum coupling in wave functions and 
E2 transition matrix elements. However, they are severely restricted by taking into account only a 
small part of many-particle many-hole excitations (single-j orbits and s,fi-bosons, respectively). The 
present approach includes all the excitations in the single-particle orbits in deformed nuclei within a 
given excitation energy. 

After describing the formulation in SectJ|, we present numerical results for ^®®Yb in Sect.^. 
We investigate in detail the excited states in the warm region of the nuclear spectrum with very 
high spins as well as associated E2 transitions, which characterize the collective rotational motion. 
In particular, we look into the properties of rotational damping and compare them with available 
experimental data. In Sectj4[ we study in detail high multipole components of the two-body effective 
force and compare the surface delta interaction with the pairing plus quadrupole-quadrupole force. 


2 Formulation 

2.1 Single-particle basis 

We consider deformed nuclei with stable prolate shape. In such a nucleus, rotational bands observed 
by the discrete gamma-spectroscopy are often interpreted in terms of the single-particle excitations 
in cranked mean-field models where the collective rotational excitation is described in a semi-classical 
way while the intrinsic particle excitations are described quantum mechanically. For example, the 
cranked Nilsson model is widely used and gives a realistic description of the single-particle excitation 
modes in the high spin yrast region MM- Let us assume that the potential surface at a given spin 
displays a deep minimum in the space of the deformation parameters for a prolate deformation. In 
such a case, it should be possible to have intrinsic excitations of many particles and holes (or many 
quasiparticles) on top of the stable deformed mean field state. The excited states above the yrast 
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line are to be built out of these excited configurations. 

In keeping with this picture, we formulate a shell model in which the single-particle basis is repre¬ 
sented by the cranked Nilsson model and the many-particle many-hole (np-nh) excitations associated 
with the single-particle potential are taken into account to represent the intrinsic excitations in ro¬ 
tating nuclei. The collective E2 transitions are assumed to keep the intrinsic structure, that is, each 
of the np-nh configurations forms a rotational band. Since the level density of the np-nh states is 
high, the residual two-body force causes mixing among the np-nh excitations, that is, mixing among 
the rotational bands. The configuration mixing can be taken into account explicitly by solving a 
shell model Hamiltonian which combines the cranked Nilsson mean-field and the residual two-body 
force 

In the following, we describe the details of the model. Let us start with the single-particle basis. 
The single-particle Hamiltonian in the present theory is given by the cranked Nilsson model 

hcrank — hNilsson ( 1 ) 

in which the protons and neutrons move in a deformed Nilsson potential which rotates around the 
axis with the largest moment of inertia (taken conventionally as the x axis), with uniform rotational 
frequency uj. We shall not include the pairing potential in the mean-field Hamiltonian. Note, how¬ 
ever, that pairing correlations may be partly taken into account through the configuration mixing 
caused by the residual interaction. Moreover, we shall be mainly interested in levels at very high spins 
I ^ 30, where the pairing gap is expected either to be small or to vanish, due to the Coriolis anti¬ 
pairing effect 11^ . Inclusion of the pairing potential in hcrank may improve the description 

of near-yrast rotational bands at low spins. However, use of the quasiparticle scheme may lead to 
incorrect description of highly excited states since many-quasiparticle configurations accompany spu¬ 
rious components related to the nucleon number violation. We use the Nilsson potential |2^ with a 
single-stretched P term and an Is force whose parameters are given in fl^ . The quadrupole and hex- 
adecapole deformation parameters of the Nilsson potential (e, € 4 ) should in principle be determined 
selfconsistently by minimizing the potential energy surface. The numerical calculations presented in 
the present paper were carried out with the deformation parameters (e, 64 ) = (0.255, 0.014) for ^®®Yb 
, according to Ref, . The single-particle basis keeps the signature and parity quantum numbers. 

The eigenstates of the cranked Nilsson Hamiltonian, Eq.(|^, define the adiabatic single-particle 
orbits. The eigenrouthians and the single-particle wave functions vary with the rotational frequency 
Lo. From time to time pairs of orbits cross each other as a function of ui (See fig.^. If we magnify the 
crossing among the orbits with the same quantum numbers, it is seen that the two orbits “repel” each 
other. The nature of the two orbits interchanges across the crossing. This causes an abrupt change 
in the single-particle orbits as a function of the rotational frequency. This abrupt change makes it 
difficult to define the rotational bands on the basis of the single-particle orbits given adiabatically. 
To remedy this problem and guarantee the smooth change as a function of w, a diabatic basis for 
the single-particle orbits is constructed in place of the adiabatic eigen solutions. The method 
of constructing the diabatic basis is described in the Appendix A. The resulting diabatic routhian 
spectrum {e'(a;)} is also shown in Fig|l[ In the following, we adopt the diabatic single-particle basis 
with routhians {e'(a;)}, single-particle wave functions {ipi{uj)} and the associated angular momentum 
expectation value 

2.2 Unperturbed rotational bands 

Once the single-particle basis is defined, the shell model basis of the many-body system is obtained by 
filling N neutrons and Z protons in the diabatic cranked Nilsson orbits. Based on the configuration 
where all the orbits up to the Fermi surface are occupied, one can generate many-particle many-hole 
(np-nh) excitations, which form shell model basis states for the excited states above the yrast line. 
The basis configurations, labeled by at spin I [| are given by 

w))= n (2) 

occupied i in fi 

^In this paper, the spin is measured in units of h 
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Figure 1: The cranked Nilsson single-particle routhian spectra for (a) neutrons and (b) protons 
with deformation parameter (e,€ 4 ) = (0.255,0.014). The different kinds of curves denote different 
parity and signature; (tt, cc) = (-I-, 1/2), (-I-, — 1/2), (—, 1/2) and (—,—1/2) for solid, dotted, dashed 
and dot-dashed curves, respectively. The circles indicate the positions where the pair-wise repulsions 
take place. Both adiabatic and diabatic basis are plotted while they mostly overlay with each other. 
In the right panel, we show a magnified portion of the neutron routhian spectrum so that difference 
between the two basis becomes visible. 
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Figure 2: The unperturbed rotational bands in ^®®Yb with (tt, a) = (+, 0). The lowest one hundred 
levels are plotted for each spin with small horizontal bars. A reference energy /(/ + 1)/2J with 
J = 76 MeV“^ is subtracted. 


where a| represents the creation operator for the diabatic cranked Nilsson single-particle wave func¬ 
tion ipi{u!i) occupied in this configuration. The rotational frequency ujj corresponding to the spin 
value I is calculated by imposing the condition 

{J:,){UJI)=I (3) 

on the average value of the projection (J^:) of the angular momentum on the rotation axis. Here this 
quantity is calculated by using thermal averaging as 

(Jx) = : ( 4 ) 

i 

where jx,i{i^) is the expectation value of jx of an orbit i, and = {1 -I- exp (e'(a;) — X)/T} is the 
thermal occupation probability. The temperature parameter is chosen as T = 0.4 MeV so that it 
corresponds to thermal excitation energy U ^ 2 MeV relevant to the excited states under discussion. 

The energy of the basis states in the laboratory frame is given by the Strutinsky renormalization 
method according to the prescription commonly used in the cranked Nilsson model |^, 

Namely, the energy A^(/) of the basis state at spin I is given by 

Ef,{I) = (/) - (I)+E^^{I) (5) 

Here the shell correction energy — Esmooth^-^ jg calculated microscopically with use of the 

diabatic single-particle basis. The bare energy E^'-^^{I) is defined by £’^***(7) = E'^{lo) luJx^^{lo) 
with a constraint Jx,^{to) = 7 on the rotational frequency u. Here E'^{lo) and Jx,p,{uj) are the total 
routhian and the expectation value of the angular momentum Jx-, respectively. Utilizing the diabatic 
single-particle routhian basis, E'^{uj) and Jx,p,{uj) change smoothly as a function of the rotational 
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frequency, and it is possible to approximate the bare energy in the laboratory frame as 





= e'(w/), 

occupied i in 

(7) 


— ^ ' jx,ii^l)j 

occupied i in 

(8) 

J?' 

du) 

5 

UJI 

(9) 


by means of the second order extrapolation from the reference rotational frequency loj |p^ ^ . The 
deviation \ Jx,ti{^i) — I\ in the angular momentum expectation value is less than 5 at spin / = 50 for 
most configurations in the present calculation. The Strutinsky smoothed energy is given by 

= ( 10 ) 

i 

and 

= I, (11) 

i 

with the smoothed occupation number fii Q. The rigid-body rotational energy = /(/ -|- 

l)/2Jrigid is calculated for the given shape of the potential {Jrigid = 80.3 MeV“^ for ^®®Yb ). Since 
the single-particle orbits have good parity tt and signature a, the shell model basis states also conserve 
the same quantum numbers. 

The energies {if^(/)} and the many-body wave functions {l/r)/))} define the shell model basis 
at a given spin / and parity tt. It is assumed that to each np-nh configuration there corresponds a 
rotational band. It is possible to make this correspondence uniquely and to follow the evolution of 
the band as a function of / in a continuous way because the basis states depend smoothly on the 
rotational frequency due to the diabatic construction of the single-particle basis. Such rotational 
bands are shown in Fig. I for i68Yb . 


2.3 Residual two-body interaction and shell model diagonalization 


Since the mean-field potential is represented by the Nilsson deformed potential, it is the residual part 
of the two-body nuclear effective force that is to be taken into account as the shell model two-body 
interaction. In order to separate the residual part, we utilize the fact that, given a shell model 
configuration (a determinantal many-body state), any two-body force is decomposed unambiguously 
into the mean-field and the residual parts. The residual two-body interaction V).es associated with 
the reference shell model configuration, denoted by /J,ref, is given by subtracting the mean-field part 
I4if from the two-body interaction; 


b-es = V2body — bnf + (V2body) 

Here the two-body interaction V 2 body is expressed as 

b2body ~ ^ 

ijkl 


( 12 ) 

(13) 


with Vijki = — 'ijjj{l)'ipi{2)\v{l,2)\'ijjk{l)i’i{2)) being the anti-symmetrized matrix ele¬ 

ments of the two-body force for the single-particle orbits i,j, k, and 1. The mean-field part bnf and 
a constant expectation value (V^body) is given by 


bnf 


Vij — Vikjk , 

(14) 



occupied k in /^re/ 


(V2body) 

= 5 i: 

kj k 

(15) 


occupied jk in 

f-^ref 
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In this way, we could construct the residual interaction so that it has no diagonal matrix elements 
for the reference configuration, i.e. (^ref iKesl^ref) = 0. In the case of a even-even nucleus, we select 
as the reference /Zref the configuration in which the single-particle routhian orbitals are filled up to 
the Fermi level at low rotational frequency. We then take the diabatic continuation of the same 
configuration to extend to higher rotational frequencies. The residual interaction thus determined is 
essentially independent on spin or rotational frequency. 

As the effective two-body effective force, we adopt the surface delta interaction (SDI) with 
standard strength Vq = 27.5(A |Q(See the Appendix B for details). The SDI has been used widely 
in the shell model descriptions of nuclei in a wide mass range from light nuclei to deformed rare-earth. 
The delta force acting only at the nuclear surface displays essential features needed for describing 
low-lying excitations. Note that it contains the the pairing force in the particle-particle channel 
as well as the multipole-multipole forces in the particle-hole channel. In a later section, we will 
investigate this two-body force by changing its force strength as well as by comparing it with other 
effective forces. 

Given the residual two-body interaction, the shell model Hamiltonian for the basis states {|^(/))} 
at spin I is constructed as 

V + (16) 

where = (/r(/)|Vres|M^(-f)) are the matrix elements of the residual two-body interaction. 

The energy eigenstates of the shell model Hamiltonian are obtained by numerical diagonalization. 
The resultant eigenstates are admixture of the shell model basis states, that is, the unperturbed 
rotational bands ^7s, 

|a(/))=^X;(/)|Ai(/)) . (17) 

In carrying out the diagonalization, truncation of the basis states is necessary. For that purpose we 
first truncate the single-particle routhian orbitals so that only those located within 5Emax above and 
below the Fermi surfaces are treated as active orbits. In the numerical calculations presented below, 
we choose SEmax = 3 MeV, which corresponds to 20-30 orbitals for protons and for neutrons. We then 
consider all the np-nh configurations whose routhian excitation energy 5E' = '^particles ^i~'l2hoies 
lies within the interval SEmax- This generates about 2000-3000 np-nh configurations for each most 
of which have Iplh to 4p4h character while there is only a little contribution of 5p5h under the energy 
cutoff. The dimension of the basis is still too large to carry out systematic numerical diagonalization. 
To avoid this difficulty, we calculate the energy {Ep{T)} of the basis states {|/.t(^))} as well as the 
diagonal contribution to the basis energy of the residual interaction V for each configuration. 
The basis configurations are sorted up according to the energy Ep{I)+V containing the diagonal 
residual interaction, and the lowest 10^ basis states are included for the shell model diagonalization. 
The diagonalization is done separately for each . 


2.4 E2 transition strengths 


Using the assumption I I^, the electromagnetic transition matrix elements can be calculated 
within the cranking model In order to calculate the stretched E2 transitions, we neglect minor 
contributions of the non-collective E2 transitions so that only the E2 matrix elements connecting 
the same intrinsic configuration are taken into account. Assuming that the nuclear shape is stable 
and independent of the configuration, the E2 operator then becomes proportional to a constant 
quadrupole moment Qo] 




{^l{I)\MiE2,X = 2,^i, 


15 

1287 r 


QoS, 


mm' > 


2)|/i'(J-2)) 


(18) 

(19) 


where Spp> = 1 if the configuration |^(/)) is the diabatic continuation of \^'{I — 
for other configurations. 

The stretched E2 transition probability among the energy eigen states a at / 
calculated as 

B{E2, al ^ a'/ - 2) = 


2 )), and S^^^ = 0 
and a' at / — 2 is 
( 20 ) 
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with amplitude 


(21) 


M 

Hereafter, we often use the normalized E2 strength 

Sal,a'1-2 = -^a/,a'/-2 ) (22) 

which satisfies Sai,a'i -2 = 1, for E2 decays from the level a at I. The normalized strength for 
a given transition coincides with the normalized transition probability or the branching ratio for the 
decay from a at / to a' at / — 2, neglecting the factor. 

In the experiments, only strong transitions are observed as discrete peaks in the gamma-ray 
spectra and the rest of the transitions shows up as quasi-continuum spectra which contain transitions 
summed over many states. Eor such situation, it is useful to represent the E2 transition properties by 
means of the distribution function of the strength, or the strength function. The strength function 
for the stretched E2 decay from the levels at I to the levels at / — 2 is given by 

Si{E^) = Y, SaI,a'I-2fal5{E^ - Eal + Ea'I- 2 ). (23) 

olcx! 

where ]ai is the feeding probability of the level a at 1. It becomes sometimes useful to define a 
strength function 

S\,a{E,~f) = 'Yj SaI,a'I-2S{Ej — Eal + Ea'I- 2 ), (24) 

a' 

for stretched E2 gamma-decays from a given specific level a at I. Average of Si,a weighted with the 
feeding probability is equivalent to the strength function defined by Eq.(^). 


3 Results for 

3.1 Mixing of t^p-nh configurations 

The calculations presented below were carried out for the deformed rare-earth nucleus ^®®Yb . For this 
nucleus there exist experimental data from the analysis of quasi-continuum gamma-spectra as well as 
data from discrete-peak spectroscopy identifying the rotational bands up to around / ~ 40 |3^, |3^ . 
The potential energy surface in this nucleus has a stable minimum at prolate deformation e ~ 0.25 
up to / ~ 60 @11,13]. and many of the observed rotational bands are described as independent 
particle excitations in the cranked mean-field ||3^, |l^. This validates the basis assumption of the 
present model. We have also made calculations for other nuclei near ^®®Yb obtaining similar results. 

The shell model diagonalization is done separately at each for / = 20 — 61. The energies of the 
levels calculated within the lowest 10^ basis states cover an interval of ~ 3 MeV. The level density of 
these states is plotted in Figj^ for a few spins as a function of the excitation energy U measured from 
the lowest state for a given . It is compared with the Fermi gas level density with fixed signature 
and parity, appropriate for the cranked mean field |@], 


Pfg{U) = 


*U * exp 2\/aU 


(25) 


where a is the level density parameter. It is seen that the calculated level density increases expo¬ 
nentially as a function of the excitation energy, and follows approximately the Fermi gas formula up 
to about 2.5 MeV above the lowest state. The level density parameter fitting the calculated level 
density is about a ~ A/10 MeV“^. This agrees with the standard estimate based on the average 
single-particle level density in the harmonic oscillator model. The decrease of the level density at 
U ^ 2.5 MeV arises from the truncation of the basis states. 

As a consequence of the configuration mixing caused by the residual interaction, the energy 
eigenstates {|a)} are admixture of many unperturbed states. In other words, the unperturbed states 
are spread over the energy eigenstates. The mixing is strongly dependent on the internal excitation 


8 




Figure 3: The level density for the calculated energy eigenstates with I'^ = 30'*", 40'*", 50'*' in 
is plotted in the left panel as a function of the excitation energy measured from the lowest state 
with same . The Fermi gas level density (long dashed lines) is also plotted with the level density 
parameter a = A/8,^/lO, A/12 MeV“^. The thick lines cover the energy region containing the 
lowest 300 levels while the short dashed lines continue above. In the right panel, the mixing number 
Timixing Is plotted for the calculated I'^ = 50+ states. The lowest 300 states are labeled with a thick 
asterisk. 


energy U of the states since the level density increases exponentially with U. A measure of the 
configuration mixing can be defined by means of a quantity 

n„.ixing(a/)= > (26) 

which counts effectively the number of basis states s participating to form an energy eigenstate 
|a) = Ia^)' l^ote that, when all the states /i have the same probability |X“| , rimixing foves 

the exact number of admixed states. This quantity, called the mixing number, is plotted in Fig.^^b) 
for a typical case at = 50+ as a function of the excitation energy U measured from the yrast 
state. For the states near the yrast (17^1 MeV), the mixing number is less than about 3, which 
indicates that each of these states is composed of a few dominant configurations. As U increases, the 
mixing number increases steeply. For the states at C/ ~ 2 MeV, several tens of np-nh configurations 
contribute to form an energy eigenstate having rather complex wave functions. It is also to be 
noticed that the mixing number fluctuates significantly state by state. There are several states with 
low mixing number rimixing < 3 even at U = 1.5 — 2 MeV where most of the others are strongly 
mixed with rimixing ^10. 

The effect of the truncation is visible both in Fig.^(a) and (b) for the levels at 17 ^ 2.5 MeV. We 
find that the lowest 300 levels located at U ^ 2.3 MeV are stable against the truncation as far as the 
level density and the gross features of the mixing number are concerned. 

3.2 Onset of rotational damping 

The calculated energy levels are plotted in Fig.^ with little horizontal bars for (+, 0) states (positive 
parity and even spin) in ^®®Yb . In this picture, a reference rotational energy is subtracted from 
the calculated energy so that the vertical axis represents roughly the internal excitation energy 
of the nucleus. The solid lines connecting the energy levels represents strong E2 transitions. We 
use the convention that the E2 decay is “strong” if the associated normalized strength satisfies 
S > l/-\/2 = 0.707. Weaker transitions with 0.5 < S' < 0.707 are displayed with dashed lines. It is 
seen that the strong E2 transitions (solid lines) form sequences of levels which are aligned regularly 
along a parabola like curve. Such sequences of levels represent rotational band structures. Most of 
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Figure 4: The calculated levels in with (+,0) are shown with small horizontal bars. A 

reference energy /(/+ 1)/2J with J = 76 MeV“^ is subtracted. The stretched E2 transitions which 
have the normalized strength Sai^a'i -2 larger than 0.707 are plotted with solid lines connecting 
initial and final levels of the transitions. Transitions with normalized strength between 0.5 and 0.707 
are plotted with dashed lines. 


the rotational bands lie in the region near the yrast, while the levels at higher excitations energy as 
a rule do not form band structures. 

To indicate the E2 transition properties more precisely, we calculate distribution of the E2 tran¬ 
sition strengths Eq.(^J), defined for gamma-decays from individual levels a at . They 

are shown in FigJ^ for the lowest levels with I'^ = SO'*'. The E2 strength associated with the first 
30+ level (which is the second lowest at spin 30) is concentrated in a single component feeding the 
28+ level with strength exhausting more than 95% of the total strength. The 30J, 30;^ and 30J 
levels show essentially the same E2 distribution except slight difference in the gamma-ray energy 
for the dominant transition. The level 30^, with an internal excitation energy of about 700keV, 
displays a completely different E2 strength distribution, being fragmented over several transitions, 
each of which carries a rather weak strength. The E2 strength associated with the decay from the 
30^, 30^, 30^ and 30^ levels shows a similar fragmentation. The fragmentation of the E2 strength 
increases as the excitation energy increases. Figure ^ displays the quantity Si^a{Ej) associated with 
the levels 30^3 and 30^4 lying at f/ ~ 1.5 MeV and for the levels 30+gQ and 30+g4 lying at ^ 2.0 
MeV. At ~ 1.5 MeV, the E2 strength distribution has about ten branches, while the number of 
branches becomes much larger at t/ ^ 2.0 MeV. The E2 strength is distributed within the range 
E^ ~ 0.7 — 1.0 MeV, with centroid at ~ 0.85 MeV and with a width of about 150 keV. 

The fragmentation of the E2 strength function is the rotational damping phenomenon Q . Aver¬ 
aging the E2 strength distributions over many states (Figj^) produces smooth profile for the strength 
distribution. The centroid of the distribution increases with the spin as expected from the relation 
{E^) ~ 2TiuJrot ~ 2// J, J being the average moment of inertia. The width of the distribution is about 
100 keV for I ^ 30, and 250 keV for I ^ 50. This approximately agrees with the estimate for the 
rotational damping width F^ot = 4Atu = 0.14(//40)C/^/"‘ MeV in Ref.[|[. To be more precise, how- 
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Figure 5: The distribution Si^a{E-y) of the stretched collective E2 decays from the lowest 9 levels 
with = 30+. The branching number Ubranch defined by Eq.(p^ for the E2 decay and the relative 
excitation energy U (in unit of MeV) are put for each level. The bins for the transition gamma 
energy Ej have a width of 6 keV. 


ever, the strength distributions are not necessarily represented by a simple Lorentzian or Gaussian, 
but exhibit structures which reflect specific alignment properties of the underlying single-particle or¬ 
bitals. In this particular nucleus, there exist two significantly aligned proton orbitals near the Fermi 
surface, 7r/ig/2 and 7 ^ 113/2 (see Fig||). These aligned proton orbitals cross the Z = 70 Fermi surface 
at w ~ 0.4 MeV or / ~ 40. The peak at ~ 1.0 MeV of the / = 40 ^ 38 strength distribution 
consists mainly of the transitions involving one of the aligned 7r/ig/2 and 71113/2 orbitals while the 
little bump at E^ ^1.2 MeV arises from the components containing neither Tr/ig/g nor 71113 / 2 . An 
effect of the aligned proton orbits is also seen in Fig.^ as changes in slopes of the rotational band 
structures near the yrast at / ~ 40 (see also the later discussion). 

It may be useful to define the onset energy where the rotational damping sets in. In order to 
quantify the onset of rotational damping, let us utilize the branching number il, EM 


77'branch(o^) 



(27) 


which counts effectively the number of the E2 branches for decays from a level a at / to levels 
at 1 — 2. This quantity is analogous to the mixing number introduced in Eq.(p^. Figures ^ and 
H also show the branching number for each level. As seen in Fig.||, the onset of damping can be 
characterized by a condition Ubranch > 2 implying an E2 decay with more than 2 branches. The 
dependence of the branching number with excitation energy is depicted in Fig.|[ The branching 
number increases exponentially with internal excitation energy. With the criterion Ubranch > 2 for 
the onset of rotational damping, the onset energy is read from Fig.^ to be 17 ^ 800 keV above yrast. 

It should be emphasized that, although the onset energy thus defined tells approximately where 
the rotational damping sets in, the transition from the region of rotational bands to the region of 
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Figure 6: The strength distribution Si^a{Ej) for the stretched collective E2 decays from the 53-rd, 
54-th, 180-th, and 181-st excited levels with I'^ = SO’*'. 



(MeV) 


Figure 7: The strength function Si{E^) for the E2 strength for the decays / = 30 ^ 28, / = 40 ^ 38 
and / = 50 ^ 48, summed over the lowest 200 states for both parities with a constant feeding 
probability /«/. 
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Figure 8: The E2 branching number Ubranch for the = 30+ (crosses), 50+ states (circles) as a 
function of the excitation energy U of the states measured from the lowest one. The histogram shows 
ribranch averaged for energy bins. 


rotational damping does not take place sharply at the onset energy, but rather develops gradually 
as the excitation energy increases. In fact, Fig.^ shows presence of many short rotational hands for 
which the strong E2 transitions (solid or dashed lines) continue for only a few to several steps. They 
lie from the border region ?7 ^ 1 MeV of the onset of rotational damping up to a region of much 
higher level density with ?7 ^ 1.5 MeV. Here the short band structures are surrounded by levels 
which do not have any strong transitions. This indicates that the rotational band structures partly 
remain even in the region of the rotational damping. Presence of such scars of rotational bands 
|^l| is also displayed in Fig.^ from which it is seen that there exist levels whose branching number 
ribranch IS Smaller than 2 or 3 even at high excitation energy U ~ 1.5 MeV. It is noticeable that 
the E2 strength distribution in the transition region is fluctuating quite irregularly from state to 
state as the FigsJ^ and || indicates. Two panels in the left (right) hand side show quite different 
fine structures, even if they are produced from neighboring levels 30^3 and 30^4 (SO+gQ and SO^^g^), 
corresponding to excitation energies which are the same within just 8 (2) keV. 

3.3 Number of rotational bands 

Since the rotational band structures do not survive in the region of high internal excitation energy, 
there exists only a finite number of discrete rotational bands in a single nucleus. This feature can be 
utilized to study experimentally the onset of rotational damping. Through the fluctuation analysis 
method i0 it is possible to extract an effective number of gamma decay paths from the double¬ 
coincident X spectrum. When this method is applied to the ridge structures in the spectrum, 
the effective number of paths essentially corresponds to the number of discrete rotational bands. 

Let us define a corresponding quantity within the framework of the present calculation to be able 
to compare our results with the experimental findings. A straightforward definition is just to count 
the number of strong E2 transitions at a given spin with strength S larger than a given threshold 
Sthr, for which one may use Sthr = l/v^ = 0.707 to define levels forming discrete rotational bands as 
adopted in plotting Fig.^ Keeping in mind the previous analysis for the onset of damping, another 
definition may also be used in terms of the branching number with criterion Ubranch < 2 for levels 
forming rotational bands. Note that these two criteria are not very different, since each state which 
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Figure 9: The calculated number of bands with length more than or equal to 2 is compared with 
the experimental effective number of paths extracted from the first ridge of {E^i,E^ 2 ) spectra in 
lesyb . The solid line is calculated with the criterion Ubranch < 2 while the dashed line is defined 
with S > 0.707. The horizontal axis denotes the average gamma-ray energy E^ = {E^i + E^i)l2. 


has a very strong transition with S > 0.707 necessarily has a branching number Uhranch < 2. The 
fluctuation analysis method is often applied to the first ridge located within the interval E^i — E ^2 = 
±4/J in the E^ x E^ spectra (J being the moment of inertia), which is formed by two consecutive 
E2 gamma-rays. Correspondingly we count the number of rotational bands which satisfy the above 
criteria at least over two steps / -|-2^/^/ — 2ofE2 decays. In practice, we count all the states 
which satisfy ntranch < 2 or S > 0.707 for both decaying transitions and feeding ones, and sum 
the numbers from four sets of parity and signature = Iq , Iq , {Iq + 1)^, {Iq + 1 )~ for a given 
representative spin Iq. The calculated number of bands can be plotted as a function of spin Iq or 
the average gamma ray energy Ej = {E^i -|- E^2)l2. Figure || shows the calculated number of bands 
using the two criteria and compares them with the experimental effective number of paths associated 
with the first ridge. The two criteria give essentially the same number of bands around 30. It is 
noticed that the theoretical calculations and the data agree quite well in all the y-ray energy range 
E^ ~ 0.8 — 1.2 MeV corresponding to spins J ~ 28 — 48. 

More accurate comparison between the theory and the experiments requires evaluation of the 
feeding probabilities. This is possible by making a simulation of the whole gamma decay cascades 
combining microscopically calculated levels and E2 transitions with statistical description of El 
decays since the present model can describe levels up to about 2 MeV above yrast line at high spins, 
where most of the gamma decay cascades of a warm rotating compound nucleus proceed. Such 
a microscopic simulation has been developed recently]^. From the E^ x E-^ spectrum produced 
from the simulation, the effective number of paths is extracted in the same way as the experimental 
analysis. In this manner, we can make comparison which does not depend on the particular choice 
of the threshold for Stur or nh^anch used for the definition of the number of bands. It is found that 
the effective numbers of paths both in the simulation and in the experiments for ^®®Yb agree quite 
well 1 ^ . 

As discussed above, the calculations predict the presence of scars of rotational bands, that is, the 
“short” rotational band structures connected by strong E2 transitions for only a few decay steps. It 
may be possible to check this feature in the experiments by looking into the second and higher order 
ridges in the E^ x E^ spectra. As the two consecutive gamma-rays from a rotational band fall on 
the first ridge at E-^i — E ^2 = 4/J, two gamma-rays of the first and last decay steps in n consecutive 
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(MeV) 

Figure 10: The calculated number of bands with length more than or equal to 2, 3, 4 and 5 steps, 
plotted with solid, dotted, and dot-dashed lines, respectively. In this figure, the criterion S > 0.707 
is used to define the rotational bands. The horizontal axis denotes the average gamma-ray energy 
E^ = {E^i+E^2)/2. 


transitions should fall on the (n — l)-th ridge at E^i — Ej 2 = 4{n — 1)/J. Thus, for example, 
the effective number of path extracted from the second ridge tells the number of rotational band 
structures which continue for at least three consecutive steps. Since the rotational bands continuing 
only for two steps do not contribute to the second ridge, we may expect the effective number of paths 
for the second ridge to be smaller than that for the first ridge. This is illustrated in Figure ^ In 
this figure, we count the number of rotational bands defined by the criterion S > 0.707 for various 
lengths of decay steps. In plotting the number of rotational bands, we used the average gamma-ray 
energy for the two gamma rays, 1 + 2^1 and I — 2n + 2^I — 2n, for the rotational bands with 
> n decay steps. The calculation indicates that a considerable part of the first ridge comes from the 
scars of rotational bands at ” high” heat energy. It also suggests that the effective number of path in 
the higher order ridges may be significantly smaller than in the first ridge. In fact, the experimental 
data for the higher order ridges also indicate this tendency [l^ . In order to make a quantitative 
comparison, however, a more careful analysis using e.g., the microscopic simulation of Ref. |^^ is 
required. 

3.4 correlation 

Since the quasi-continuum E^ x E^ spectrum displays the characteristic ridge-valley structure, the 
shape of the spectrum provides us with important information about the rotational damping i0- 
The quasi-continuum analysis often deals with projection of the Ej x E^ spectra upon the E^i — Ej 2 
axis. A quantity which is related to the projected spectra is the strength distribution function of 
the coincident two E2 transitions. In particular, the two-gamma strength function for consecutive 
E2 decays I + 2^I^I — 2is important since it characterize the shape of the first ridge as well as 
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Figure 11: The two-step strength distribution S 2 {Eyi — £^ 2 ) of consecutive E2 transitions for 
different bins of excited levels, defined for the lowest 10 levels, the 11 -st to 20 -th levels, the 21 -st to 
40-th, and 41-st to 60-th for each I'^. The strength is averaged for the spins / = 30 — 51 in order 
to get enough statistics to produce smooth prohle. The strength distribution is also subdivided 
with respect to the product strength pip 2 = Sai+ 2 ,a'iSa'i,a"i -2 of individual transitions. The 
approximate excitation energy of the bin measured from yrast is put for each bin for the sake of 
reference. 


of the central valley. The two-gamma strength function is given by 

S2{E^I ^^ 2 ) — ^ ^ 1 Sa'JJ— 2 /a'E^2 4 - E^ ^y /— 


where 

E^^a. —— E(yJJ^2 E(yf j 

Ej^a' ^a" EyyfJ E(yff J — 2 

and the fai is the feeding probability of the level a at spin /. 


(29) 

(30) 


The calculated two-gamma strength function S 2 {E^i — E-y 2 ) is depicted in Figure |lll Here an 


average is taken over the levels within a spin interval / = 30 — 51. In order to study the excitation 
energy dependence, the sum over a' is divided into bins including 10 or 20 levels for each I'". The 
feeding probability fa'i is put equal for all the levels. The lowest bin (^1—^10) covers approximately 
the energy region U = 0 — 0.9 MeV, thus it mostly contains the E2 transitions associated with the 
rotational band structures near the yrast line. The higher bins cover the region above the onset of 
rotational damping. 

For the lowest bin, the most characteristic feature is the presence of a sharp peak located at 
Eji — E ^2 ~ 70 keV with a width of about 30 keV, which corresponds to the first ridge observed in 
the E^ X E^ spectrum. In order to show the contents of the peak, the two-gamma strength function 
S2{E-yi — E~f2) is subdivided by putting cuts on the product strength Sai+2,a'iSa'i,a"i-2 = P 1 P 2 of 
each transition (a, J-|- 2) ^ {(3,1 -\-2) —> ( 7 , / — 2). The sharp ridge in the lowest bin consists mostly 
of the strong transitions satisfying pip 2 > 0.5. This indicates that the sharp ridge is formed by the 
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strong transitions (p > l/'s/S) associated with the rotational band structures. The peak position 
of the sharp ridge is related to an average dynamic moment of inertia J of these bands through 
E^i — E ^2 = 4/J, and the peak width originates from fluctuations in the moment of inertia among 
different bands. 

For the higher bins (#11 — #20, #21 — #40, #41 — #60), the profile of the two-gamma strength 
function is quite different from that of the lowest bin in many respects. Note that there is only a 
little contribution of strong E2 transitions (piP 2 > 0.5) associated with the rotational band structure. 
This is because the rotational damping causes fragmentation of the E2 strengths with a spread in 
gamma ray energy of order of Tj-ot. Moreover, the spectrum shows a two component profile with 
wide and narrow widths of about 300 and 80 keV. The intensity of the wide component increases 
with the average energy of the bins. This component is easily understood as a consequence of 
the rotational damping. If two consecutive transitions were uncorrelated, the two-gamma strength 
function S 2 {Eji — E^ 2 ) would simply become a convolution of two single-gamma strength functions 
Si{Ej) for the consecutive steps, producing a smooth distribution with a width a factor of 2 larger 
(assuming a Lorentzian shape) than that of the single-gamma strength function Si{E^), which 
displays a width of about 100-200 keV (see Fig.^). 

The narrow component has a width of about 80 keV, whose intensity decreases with the energy of 
the bin. The origin of the narrow component is visible in Fig.^by subdividing the strength function 
with respect to the product strength pip 2 of individual transitions. It is noticed that the transitions 
building up the narrow component have relatively large strength compared to those composing the 
wide components. Taking the third bin (#21 — #40) as an example, the narrow component consists 
of transitions with 0.03 < piP 2 < 0.3 while weaker transitions piP 2 < 0.03 contribute only to the 
wide component. These transitions in the narrow component are much stronger than the average 
strength {piP 2 ~ iXrotp)~^ ~ 0.01) expected for damped transitions. 

This analysis indicates that the correlation E^i — E ^2 ~ 4/J associated with the rotational band 
structure still remains in the consecutive E2 transitions even after the rotational damping sets in 
and that the narrow component has the same origin as the “scars” of the rotational bands discussed 
above. This suggests that the width of the narrow component is related to the spreading width 
which represents the extent of the configuration mixing of the rotational band structures. 

It is also noticed that the narrow component is intense even at the fourth bin (#41 — #60) 
corresponding to excitation energy 17 ~ 1.5 MeV above yrast. It becomes less intense as the excitation 
energy increases, and almost disappears at C/^2 MeV as shown in Fig.^. The development of 
rotational damping is quite gradual as a function of 17, and the transition region extends from 
U ~ 0.8 MeV to 17 ~ 2 MeV above yrast line. Interestingly, it is found from the statistical analysis 
of energy level spacing and individual E2 strength |^, ^ that the E2 strength shows the Porter- 
Thomas distribution at 17 ^ 2 MeV while it deviates significantly from the random limit for 17 ^ 2 
MeV. The deviation from the Porter-Thomas distribution and the presence of the narrow component 
are related with each other since both originate from the presence of strong transitions in the region 
of the rotational damping. 

The Ej — E^ correlation only displays a weak dependence on angular momentum, as shown in 


Figure 13, The intensity of the narrow component becomes larger for lower spins although there 


exist significant narrow components even at very high spin (7^50). The smaller the spin is, the 
smaller becomes the rotational damping width (width of the wide component). When the damping 
width becomes smaller, the branching number and strength of the fragmented E2 transitions becomes 
smaller and larger, respectively. Thus strong transitions which cause the narrow components become 
more dominant at lower spins. On the other hand, the width of the narrow component is not very 
dependent on spin. This also suggests that the narrow width may be related to the spreading width, 
which is essentially spin independent in the present model. 

In the experiments, the actual E^ x E^ spectrum is not simply represented by the two-gamma 
strength function discussed above since it is formed not only by the consecutive E2 transitions but 
also by non-consecutive E2’s as well as El transitions all weighted by the feeding probability. It 
is found, however, by means of the simulation analysis that the two-component profile causes 
significant effects on the spectral shape of the E^ x E^ spectra [pT|. 
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Figure 12: Excitation energy depen¬ 
dence of the two-gamma strength dis¬ 
tribution S' 2 (E-yi — £’ 72 ) for the energy 
bins covering the first to 50-th, 51-st to 
100 -th, 101 -st to 200 -th, and 201 -st to 
300-th lowest levels at each I'^. The 
strength is averaged over spin interval 
I = 30-51. 



Figure 13: Spin dependence of the two- 
gamma strength distribution S 2 {E^i — 
Ery 2 ) for the energy bin covering the 
lowest 100 levels at each I'^. The 
strength distribution is shown for dif¬ 
ferent spin intervals; / = 22 — 29,30 — 
39,40-49,50- 59. 


3.5 Near-yrast rotational bands 

The main aim of the present model is to describe the internally excited states and the rotational 
damping phenomena in the warm region of the spectrum at very high spin. However, it is interesting 
to study how the two-body residual interaction affects the rotational bands near the yrast line as 
compared to the standard cranked mean-field calculations which are often used to describe those 
states. 

The calculated and observed rotational bands are compared in Figure |^. The rotational bands 
below spin / ^ 20 are not reproduced well in the present model. This is because the present 
calculation does not give a satisfactory description of the pairing correlations, which are important 
for the rotational bands at low spins. Pairing correlations may also account for the slightly lower 
moments of inertia of the experimental bands in the angular momentum region around / ^ 20 to 
40. This is taken into account by using a reference moment of inertia in Fig.|^, which is chosen 
smaller by about 10% for the experimental bands than for the theoretical ones. In the higher spin 
region with / ^ 30, the model reproduces fairly well the overall features of the lowest few bands near 
the yrast line; above the crossing at / ~ 25, the (—, l)i band becomes the yrast band, leaving the 
(-|-,0)i as the second lowest, both in the theory and in the experiment. The dominant configuration 
in the (—,l)i band is neutron Iplh excitation (t/[521]l/2)^(^[642]5/2)~^ relative to the reference 
configuration (the (+,0) band marked with R in FigjlJ). The third lowest ((+), (0)) band in the 
experiments seems to correspond to the (-I-, 0)2 bands in the model, whose dominant configuration is 
the neutron excitation (z/[521]l/2)^(:/[642]5/2)“^. Furthermore, both the model and the experiment 
show the band crossings at / ^ 40. In the calculation, the crossings involve the aligned proton orbits 
7 r/ig /2 and/or (See Fig.|^). They are, however, slightly sharper than those in the observed 

bands. This discrepancy also may be attributed to the insufficient pairing correlation in the present 
calculation (Note that Ref.|^^ discusses the smoothness of the crossing in connection with the proton 
pairing). 

In order to show the effects of the interaction, we compare in Fig. 14 with the unperturbed 
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Figure 14: The rotational band structures near yrast line in . The calculated band structure 

is shown in the left-top panel in a similar way to Figj^ subtracting the reference rotational energy 
/(/ -I- 1)/152 MeV. The right-top panel shows the unperturbed rotational bands calculated without 
the residual two-body force. See text for marks. The left-bottom panel shows the observed rotational 
bands subtracting the reference energy /(/ -I- 1)/142 MeV Q. 


rotational bands which are obtained without the residual interaction. It is noticed that the (-I-, 0) 
configuration taken as the reference (marked with R) gains energy due to the correlation caused 
by the residual interaction, which includes a part of the pairing effect. On the other hand, the 
(-|-,1) unperturbed rotational band marked with a diamond is pushed up by the residual interaction. 
This band has configuration (i^[521]l/2)^(i^[642]5/2)“^ for neutrons and (7r/i9/2)^(7r[411]l/2)“^ for 
protons with respect to the reference. Since the neutron and the proton configurations have different 
spatial density distribution (oblate and prolate along the symmetry axis, respectively), the residual 
SDI force makes this configuration energetically unfavored. The residual interaction gives, at least 
in this case, an overall improvement of the description of the near-yrast rotational bands compared 
to the unperturbed rotational bands with pure independent particle configuration. 


4 Residual interactions 

4.1 Interaction strength 

Since the rotational damping is controlled predominantly by the configuration mixing caused by the 
residual interaction, it is interesting to examine the dependence on the two-body residual force. It 
is also noticed that there is some ambiguity in the strength of the SDI (the literature value ranges 
as Vb ~ 20/4 — 30/4 MeV, see Appendix B). To examine the dependence of the results on the 
SDI strength, we performed calculations with various SDI strengths (Vb = 14/4,20/4,35/4,50/4 
MeV). We also check how the results depend on the specific features of the SDI that emphasize the 
interaction at the nuclear surface and give large matrix elements for the single-particle orbits near 
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Figure 15: Number of rotational bands calculated with various SDI strengths and with the delta 
force. The definition of the number of bands is based on the criterion that ntranch < 2 for at least 2 
steps. 


the Fermi surface. For that purpose, we also perform a calculation with a volume delta force 

v(l,2) =-vsS^xi - X 2 ) (31) 

which does not have the surface effects. For the strength of the volume delta interaction, we used 
vs^nn = vs^pp = 340 fm^MeV and vs^pn = 500 fm^MeV taken from Ref.|Q, which are evaluated from 
comparison with more realistic residual interactions, independently of the SDI strength. 

Figure ^ compares results from these forces with respect to the number of rotational bands 
taken as an measure of the onset of rotational damping. The calculated number of bands is rather 
stable for reasonable values 20/A — 35/A MeV of the SDI strength (See Appendix B), in fairly good 
agreement with the experimental effective number of path. The volume delta force also gives similar 
results. 

In contrast, if the SDI strength is reduced, e.g., to Vb = 14/A MeV, the number of bands 
increases significantly and overestimates the experimental data. This is because the configuration 
mixing effects is reduced accordingly (e.g. Umixing becomes about factor four smaller than with 
Vb = 27.5/A MeV) and the weakened configuration mixing is less capable of causing the rotational 
damping. On the other hand, the number of rotational bands does not decrease very much increasing 
the SDI strength. This feature is qualitatively explained by the role of the level densities for the onset 
of damping. According to Ref.[||, the onset energy I/q scales with the interaction strength v only 
in powers Uq ~ while the number of levels (or number of rotational bands) below Uq depends 

on Uq in an exponential manner exp\/2aUo as estimated from the Fermi gas level density. It should 
be remarked, however, that the rotational band structure near yrast line is significantly deviated 
from the unperturbed rotational bands (corresponding to the conventional cranking calculations) if 
the strength exceeds more than 35/A MeV. Such large strength may not be very realistic for the 
description of low-lying rotational bands near the yrast even though the onset of rotational damping 
is described fairly well. 

4.2 Role of high multipole components 

In addition to just the strength of the residual interaction, also the spatial properties will affect 
its ability to generate configuration mixing. Earlier pH | we have shown that a pure pairing plus 
quadrupole-quadrupole (P+QQ) interaction produces transition strength distributions and level 
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Figure 16: Distribution of the off-diagonal two-body matrix elements Vijki of SDI with Vq = 27.5 jA 
MeV, binned within intervals of width 1 keV. The angular momentum chosen for the figure is 
I = 40,41. The distribution for the low-multipole part {L < 3) and for the high-multipole part 
(L > 6) of the same SDI are also shown. The latter two have essentially the same r.m.s. value. 


spacing distributions which are very different from those generated by the SDI. The P-l-QQ force 
provides a good description of only the large matrix elements which lead to collective properties of 
low-lying levels. It is to be noticed however that the low multipole forces such as P-l-QQ usually 
have strong selectivity or selection rules with respect to the Nilsson asymptotic quantum numbers. 
Contrary, the components with high multipolarity may not have such large matrix elements nor 
cause the collective effects while they may have less selectivity. 

Figure ^ shows the distribution of off-diagonal two-body matrix elements between all possible two 
particle transitions (ij kl) conserving parity and signature associated with the interaction matrix 
elements (/r|Vres|p-0 between the basis configurations {|p.)}. The selectivity of the low multipole 
part {L < 3) of the SDI manifests itself as a surplus of very weak matrix elements, and also of 
very strong matrix elements. The restriction to very large multipoles {L > 6) is seen to lead to a 
more smooth distribution. To investigate how such a selectivity affects the configuration mixing, 
we perform several calculations in which only the high or the low multipole components of the SDI 
are taken into account. Low multipole components include L < 2,3,4, 5 in Eq.(B^) with the SDI 
strength Vq = 27.5/A MeV, and high multipoles with L > 3,4, 5, 6. Another calculation has also 
been done using the P-l-QQ interaction with the monopole pairing strength Gq = 20 /A MeV ,the 
quadrupole pairing strength G 2 = 2.Axseif which is estimated from the multipole decomposition of 
the delta force [^, and the QQ interaction strength Xnn = Xpp = ^-^XseifjXnp = ^Axseif estimated 
from the selfconsistency condition both for isoscalar and isovector components. A calculation 



r.m.s. (keV) 


r.m.s. (keV) 

all L 

19.8 



L<2 

11.3 

L>3 

16.5 

L < 3 

12.3 

L > 4 

15.7 

L < 4 

16.3 

L>5 

13.4 

L<5 

16.8 

L>6 

12.3 


Table 1: The root mean square value of off-diagonal matrix elements of the residual 

two-body force for various multipole decompositions of SDI with Vq = 27.51 A MeV. The angular 
momentum chosen for the table is / = 40,41. 
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Figure 17: The energy levels and rotational band structure for (+,0) calculated with only 

(a) the low multipole including L < 3 terms of SDI in Eq.(B.3) and (b) the high multipoles L > 4. 
The strength of SDI is Vb = 27.5/71 MeV. See Fig.^ for definition of solid and dashed lines. 


with the P+QQ interaction with doubled strength is also performed. To show the average size of 
the the two-body matrix elements, the root mean square (r.m.s.) value of the off-diagonal elements 
is listed in Table |[ Note that the r.m.s. value associated with the P-l-QQ is 12.7 keV so that it 
becomes larger than the r.m.s. of the full SDI when the force strength is doubled. 

The influence of the different multipolarity on the onset of rotational damping is evident in Fig.p^, 
in which the rotational band structures are displayed in the same way as Fig.^ for the calculation 
with low multipoles L < 3 of SDI as well as for the one with complementary high multipoles L > 4. 
With the low multipole part, most of the states shown in the figure form a large number of rotational 
band structures characterized by sequences of strong E2 transitions, apparently different from the 
calculation Fig.^ with full SDI. On the contrary, the calculation with the higher multipoles of the 
SDI gives essentially the same behavior as the full SDI except details of individual levels. 

To give a more detailed analysis, the number of bands is calculated for various choices of multipoles 
and plotted in Figjl^ It is noticeable that the low multipoles with L < 3 of SDI as well as the P-l-QQ 
interaction produce several hundreds of rotational bands, which is much larger than the full multipole 
SDI and apparently disagree with the experimental data. On the other hand, the high multipoles 
with L > 4, 5 components give essentially the same number of rotational bands as the full multipole 
SDI. These results point to that the onset of rotational damping essentially originates from the high 
multipole components L ^ 4 of the two-body residual interaction and the lower multipole components 
L ^ 3 contribute little. It should be noted that the P-l-QQ with doubled strength gives still too large 
number of bands while it has larger r.m.s. value of matrix elements than the full multipole SDI. 
The low multipoles L < 3 and the high multipoles L > 6 have the same value of r.m.s. while the 
predicted onset of rotational damping is completely different. 

The origin of this difference can be traced back to the distributions of the two-body matrix 
elements shown in Fig.|l^, where the distribution of off-diagonal two-body matrix elements Vijki is 
plotted for the SDI with all the multipoles, with the low multipoles containing ft < 3, and with 
the high multipoles L > 6 (The same plot for the P-l-QQ is given in Ref. ID)- The latter two 
have about the same r.m.s. value of about 12 keV (cf. Table 1). All the three interactions have a 
distribution which is sharply peaked at Vijki = 0 as compared to the Gaussian distribution (inverse 
parabola curve) which is expected in a random limit. In particular, the distribution of the low 
multipole interaction is significantly concentrated at small values \vijki\^ a few keV, indicating 
strong selectivity for the matrix elements of the low multipole interaction while the high multipole 
part of the SDI as well as the full SDI has more uniform distribution. On the other hand, the 
value of P(vijki) for the low multipole part {L < 3) of SDI is factor 2-3 smaller than for the high 
multipole part {L > 6) or the full SDI at around the average strength (?; ~ 10 — 20 keV) of the 
matrix elements. This implies that the basis np-nh configurations {|/r)} have less chances to interact 
with each other via the low-multipole parts of SDI than via the high-multipole parts. In fact, we 
find that the mixing number rimixing, Eq. (|^) , calculated with the low multipole part (L < 3) of 
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Figure 18: The number of rotational bands for various multipole decompositions of SDI plotted as 
a function of average gamma-ray energy (See text for details). The definition of the number of bands 
is the same as in Fig.^ The number of bands calculated with the P-I-QQ force is also plotted with 
dotdashed line. The result for the P-I-QQ with doubled force strength is plotted with dashed line. 


SDI is significantly smaller than those obtained with full SDI or the high multipole part (L > 6 ). 
For example, the average value of rimixing calculated at = 50+ is around 3 even at U ~ 1.5 
MeV, and many of the states around this energy keep a rather pure configuration with rimixing ^ 2 
(cf. {rimixing) ~ 10 and ~ 8 for the full SDI and the high multipole parts L > 6, respectively). 
Accordingly, one finds that the discrete rotational bands extend to a much higher excitation energy 
when only the low-multipole parts of SDI are included. 

One may consider several reasons for the distribution of the two-body matrix elements being 
much sharper than the Gaussian. The Gaussian distribution could be obtained if the wave functions 
of single-particle orbits had no specific structures and were randomly distributed, while such a 
random situation is not satisfied for the actual cranked single-particle orbits under consideration 
since the cranked single-particle orbits keep some characteristic structures even in the presence of 
the deformation and the cranking terms in the mean-field Hamiltonian. As an example of such 
structures, one could consider angular momentum components of single-particle orbits. Many of the 
cranked single-particle orbits can be classified as deformation aligned or rotation aligned, depending 
on whether the angular momentum vector points to the deformation axis {z axis) or the rotation 
axis {x axis), respectively. The deformation aligned orbitals, which correspond to the almost straight 
lines with small slope in Figj^, generally carry large and fairly pure K values (projection of angular 
momentum along 2 ; axis). The rotation aligned orbitals corresponding to straight lines with large 
slope are dominated by the components with large total angular momentum J (e.g. 113/2 for neutrons, 
^ii/ 2 i^ 9 / 2 i and ^ 13/2 for protons). These approximate quantum numbers may give rise to strong 
selectivity in the two-body matrix elements. This effect is expected to be stronger for low multipole 
interactions than for those containing high multipole contributions. If the angular quantum numbers 
K and J of the single-particle orbits are good quantum numbers and the two-body interaction with 
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multipolarity L is considered, there is a selection rule for the K and J of the interacting two particles 
such as \Ki — K 2 \ < L and | Ji ± J 2 I < L. This selection becomes effective for L which is smaller 
than the average J of single-particle orbits near the Fermi surface. 

5 Conclusions 

We have presented results of shell model calculations for rapidly rotating warm nuclei. The model 
describes microscopically levels and stretched E2 transitions up to a few MeV above yrast at high spin 
region I ^ 20. From the numerical calculation for ^®®Yb , it is predicted that the rotational damping 
sets in around U ^ 0.8 MeV above the yrast line. The levels near the yrast line form rotational band 
structures and the number of the rotational bands is calculated to be around 30 at a given rotational 
frequency, agreeing with the experimental data extracted from the fluctuation analysis of Ej x 
spectra. 

The model predicts some novel features of the rotational damping. The onset of rotational 
damping takes place quite gradually as a function of internal excitation energy and it shows large 
fluctuations depending on individual states. The transition region extends from 0.8 MeV to 2 MeV 
above the yrast line. Even in the region U ~ 0.8 — 1.5 MeV, there remain scars of discrete rotational 
band structures, which are characterized by short sequences of strong E2 transitions. In keeping with 
these results, the corresponding strength function associated with two-fold E2 rotational transitions 
is expected to display a two-component profile with a narrow component generated by strong and 
correlated transitions surviving in the region of the rotational damping and a wide component whose 
width is related to the rotational damping width. The origin of rotational damping can be traced 
back to the high-multipole components {L ^ 4) of the residual two-body force acting among the 
unperturbed cranked shell model states. 

The overall aim of the present work has been to extend the cranked mean field to finite temper¬ 
atures, including an effective two-body force which is capable of coupling many-particle many-hole 
excitations built upon the cranked mean field. The main advantage of the model is the inclusion of 
all excitations of the independent particle motion up to around 2 MeV above yrast line, and this 
implies a realistic behavior of the level density. As a drawback, the truncation of the shell model 
space adopted in the present paper does not allow for surface vibrations and pair-correlated states 
to be generated. However, since the effective two-body force applied includes terms which generate 
such correlations, it would be an interesting future subject to extend the shell model space and to 
study whether the shell model approach is able to describe also the correlated states at high spins. 
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Appendix A. Diabatic cranked Nilsson single-particle basis 

The adiabatic basis or the eigen solutions of the cranked Nilsson Hamiltonian sometimes show abrupt 
changes in the routhian energy {ei(w)} and wave functions {^i(w)} when the rotational frequency 
is varied. This is caused by the “repulsive interaction” associated with crossings among two orbits 


24 



with the same quantum numbers. The diabatic single-particle basis can be constructed by removing 
the unwanted repulsion at each of the crossings. 

First, we need to specify the unwanted crossings which cause the abrupt changes in the adiabatic 
basis. In order to measure the abruptness, it is necessary to give a relevant scale for the changes in 
the rotational frequency. Since the rotational band structure associated with the collective rotation is 
concerned, we consider the change in the frequency corresponding to the angular momentum change 
AI = 2 of two unit. Namely, we examine the change in {e^(a;)} and {'tjji{io)} against the variation 
from ojj to a;/_|_ 2 , where w/ is the rotational frequency corresponding to a given angular momentum 
I (See, sect j2.2| ). The relation between ui and I approximately follows I ~ Joj, J being the moment 
of inertia. 

Now the unwanted crossings are specified by the following three criteria, i) The single-particle 
wave functions at uji and w /+2 should not be very different for the diabatic basis. Given an orbit i 
at uji, an overlap condition 

I (cj/) (ca/-|_2)) I ^thr (-^■^) 

is checked for orbits j at coi+ 2 - If there is no orbits satisfying this condition, we consider the 
orbit i as interacting at w/. ii) The single-particle contribution to the dynamic moment of inertia 
= —cPe^/cpLO = djx,i/duj is calculated at ujj. If its absolute magnitude exceeds a threshold 
jthry the orbit i is considered as interacting at ioj since the moment of inertia for the configurations 
involving this orbit deviates significantly from the other configurations, iii) Under the assumption 
that only two diabatic levels interact simultaneously with a constant interaction, all the features 
of the level repulsion are specified by the interaction strength Vmt, the crossing frequency ujcrs and 
another quantity 6 which is the relative energy shift of diabatic levels under the change w/ ^ W 7+2 
in the rotational frequency. A measure of the abruptness of the crossing is given by a dimensionless 
parameter Vmt/S [^. For small VmtlS, below a threshold ethr, the adiabatic basis shows an abrupt 
change near the crossing point. Here the interaction strength Vmt and the crossing frequency uJcrs are 
extracted by searching the frequency where the distance |e((a;) — e 2 (w)| between the interacting orbit 
pair is minimum. The interaction strength Vint is a half of the minimum distance. The parameter S 
is expressed in terms of the diabatic energies of the interacting orbits, i.e., 

^ ^dia,2-t-l (^dia,z (^d+2) ^dia,2-t-l (^^-t-2))} (-^■^) 

where the expression for the diabatic energy is given by (A.5). The conditions i) and iii) are essentially 
the same as those introduced by Bengtsson 

The three criteria discussed above represent similar conditions about the abruptness, and are not 
mutually independent. In practice, we first pick up candidates with use of the first two criteria, and 
then examine them in terms of the third criterion in order to finally determine the pairs of orbits for 
which the diabatic basis should be constructed. A reasonable choice of the thresholds is found to be 
nthr = 0.933, jthr = 20 MeV“i and ethr = 1-8. 

Once the unwanted crossings are specified, we then remove the interaction among the two crossing 
orbits. In doing this, we assume that the crossings are isolated and that each crossing is represented 
by a two-level model with a constant interaction strength . For the two-level model, the relation 
between the adiabatic and the diabatic basis is analytically expressed. The routhian and the wave 
function of the diabatic basis are expressed as 


•dia,l or 2 ~ 

f{e[ - e'.^)li^ - + (e'l -k e!^)!^ 

(A.3) 

V’dia,! 

= cos{9)'ijji — sin{9)'ijj2 

(A.4) 

t/'dia,2 

= sin(6>)'!/;i -1- cos{9)'tp2 

(A.5) 

cos{9) = 

i-edi,2\/2 

(A. 6 ) 

V V 61-62 J 


in terms of the eigensolutions of the cranked Nilsson single-particle Hamiltonian, which form the 
adiabatic base. This prescription is applied to the interacting orbit pairs within a frequency interval 


25 





between ujmin and ojmax, at which the distance between the interacting adiabatic orbits becomes 10 
times the interaction strength Vmt ■ There is only small difference between the diabatic and adiabatic 
basis outside of this frequency interval, where the latter should be used. 


Appendix B. Surface Delta Interaction 

The cranked Nilsson single-particle basis is expanded in the harmonic oscillator basis (pnijm{rt0t4>t) 
= {Pt4>t) in the single-stretched polar coordinate system where the single¬ 
stretch means the scaling Xi = \J wq with the anisotropic oscillator frequencies of the Nilsson 

potential. Correspondingly, the SDI can be expressed in terms of the stretched coordinates 


?;(1, 2) = -y'(5(T; - 4)(5(rt - K) 

- i?) (5(rt,2 - i?) , 

= —V - 

rt.i rt,2 


(B.l) 

(B.2) 


where (5(nt,i2) is the delta function for the angle variables. As it is usually prescribed, the radial 
matrix elements of the SDI interaction is replaced by a constant. Then the SDI interaction acting 
on the angular and spin wave function is given 


Z,{I, = -4^Vo V£M(^t,iMnM(0t.2c/>t,2) (B.3) 

LM 

The evaluation of the angular matrix elements of the spherical harmonics is straightforward 1^ . 
The matrix elements in the cranked Nilsson orbits are then calculated directly. 

The interaction strength Vq is directly related to the strength Gq of the standard monopole 
pairing force |2^, for which the value around 20/A — 30/A MeV is often used. The requirement of 
the selfconsistency between induced mean-field and density gives an estimate Vq = 20/A MeV |^ . 
From fits in light nuclei (A ^70), it is estimated that Vb ~ 25/A MeV The analysis of the 
low-lying collective vibrations and the pairing in rare-earth deformed nuclei gives Vq = 27.5 ± 6.5/A 
MeV Thus, we adopt Vq = 27.5 jA MeV as a representative strength of SDI. 


References 

[1] G.A. Leander, Phys. Rev. C25 (1982) 2780. 

[2] B.Lauritzen, T.Dpssing and R.A.Broglia, Nucl. Phys. A457(1986) 61; 

R.A. Broglia. T. Dpssing, B. Lauritzen and B.R. Mottelson, Phys. Rev. Lett. 58 (1987) 326. 

[3] J.C.Bacelar, G.B.Hagemann, B.Herskind, B.Lauritzen, A.Holm, J.C.Lisle,P.O.Tjpm, Phys. Rev. 
Lett. 55(1985) 1858. 

[4] J.E.Draper, E.L.Dines, M.A.Deleplanque, R.M.Diamond, and F.S.Stephens, Phys. Rev. Lett. 
56(1986)309. 

[5] F.S.Stephens, J.E.Draper, J.L.Egido, J.G.Bacelar, E.M.Beck, M.A.Deleplanque, and 
R.M.Diamond, Phys. Rev. Lett. 57(1986)2912; 

F.S.Stephens, J.E.Draper, J.C.Bacelar, E.M.Beck, M.A.Deleplanque, and R.M.Diamond, Phys. 
Rev. Lett. 58(1987)2186; 

F.S.Stephens, J.E.Draper, M.A.Deleplanque, R.M.Diamond, and A.O. Macchiavelli, Phys. Rev. 
Lett. 60(1988)2129. 


26 




[ 6 ] B.Herskind, A.Bracco, R.A.Broglia, T.D 0 ssing, A.Ikeda, S.Leoni, J.Lisle, M.Matsuo, and 
E.Vigezzi, Phys. Rev. Lett. 68(1992) 3008. 

[7] T. D0ssing, B. Herskind, S. Leoni, M. Matsuo, A. Bracco, R.A. Broglia, and Vigezzi, Phys. Rep. 
268 (1996) 1. 

[8] B.Herskind, T.D0ssing, S.Leoni, M.Matsuo, and E.Vigezzi, Prog. Part. Nucl. Phys. Vol.28 (Perg- 
amon 1992) p.235; 

B. Herskind, S. Leoni,T. D0ssing,G. B. Hagemann, P. Bosetti, A. Bracco, S. Frattini, E. Vigezzi 
and M. Matsuo, Acta Phys. Pol. B26 (1995) 153. 

[9] B. Herskind, T. Dpssing, D. Jerrestam, K. Schiller, S. Leoni J. Lislie, R. Chapman, F. Khazaie 
and J.N. Mo, Phys. Lett. B276 (1992) 4; 

B. Herskind, T Dpssing, S. Leoni, M. Matsuo, N. Nica,, D.C. Radford and P. Rasmussen, Nucl. 
Phys. A557 (1993) 191c. 

[10] B. Herskind, T Dpssing, P. Rasmussen, S. Leoni, A. Bracco, R.A.Broglia, P.Bosetti,E.Vigezzi, 
M.Matsuo, preprint 

[11] S. Leoni, B. Herskind, T. Dpssing, P. Rasmussen, P. Bosetti, A. Bracco, R. Broglia, S. Frattini, 
M. Matsuo, N. Nica, E. Vigezzi, A. Atag, M. Bergstrom, A. Brockstedt, H. Carlssont, P. Ek- 
stromt, F. Ingebretsen, H.J. Jensen, J. Jongman, G.B. Hagemann, R.M. Lieder, T. Lonnroth, 
A. Maj, B. Million, A. Nordlund, J. Nyberg, M. Piiparinen, H. Ryde, M. Sugawara, P.O. Tjpm, 
and A. Virtanen, Nucl. Phys. A587 (1995) 513. 

]12] T.Guhr and H.A. Weidenmiiller, Ann. Phys. 193 (1989) 489. 

[13] J.L.Egido and A. Faessler, Z.Phys. A339(1991) 115. 

[14] S.Aberg, H.Flocard,W.Nazarewicz, Ann. Rev. Nucl. Part. Sci., 40(1990)439 

[15] Z.Szymanski, Fast Nuclear Rotation (Clarendon Press, Oxford, 1983) 

[16] M.J.A. de-Voigt, J.Dudek, Z.Szymansky, Rev. Mod. Phys. 55 (1983)949 

[17] R. Bengtsson and S. Frauendorf, Nucl. Phys. A314 (1979) 27; A327 (1979) 137. 

[18] T.Bengtsson and I.Ragnarsson, Nucl. Phys. A436(1985) 14. 

[19] S.Aberg, Phys. Rev. Lett. 64(1990) 3119; 

S.Aberg, Prog. Part. Nucl. Phys. vol.28 (Pergamon 1992) p.ll. 

[20] 1. M. Green and S. A. Moszkowski, Phys. Rev. 139 (1965)B790; 

R.Arview and S.A.Moszkowski, Phys. Rev. 145 (1966)830. 

[21] M. Matsuo, T. Dpssing, E. Vigezzi and R.A. Broglia, Phys. Rev. Lett. 70 (1993) 2694. 

[22] M. Matsuo, T. Dpssing, B. Herskind, S. Frauendorf, E. Vigezzi and R.A. Broglia, Nucl. Phys. 
A557 (1993) 211c. 

[23] A.T.Kruppa, K.F.Pal, N.Rowley, Phys. Rev. 052(1995)1818. 

[24] T.Mizusaki, T.Otsuka, and P.von Brentano, Nucl. Phys. A598(1996) 47. 

[25] J.D.Garrett, Nuclear Structure 1985 (Elsevier Science, 1985) p.lll, ed. R.Broglia, 
G.B.Gagemann and B.Herskind 

]26] Y.R.Shimizu, J.D.Garrett, R.A.Broglia, M.Gallardo and E.Vigezzi, Rev. Mod. Phys. 
61(1989)131. 

[27] Y.R.Shimizu, Nucl.Phys. A520(1990) 477c. 


27 



[28] S.G. Nilsson, C.F. Tsang, A. Sobiczewski, Z. Szymanski, S.Wycech, C. Gustafson, 1. Lamm, P. 
Moller, and B. Nilsson, Nucl. Phys. A131(1969) 1. 

[29] R.Bengtsson, S.Frauendorf and F.-R.May, Atomic Data and Nuclear Data Tables, 35(1986) 15. 

[30] T.Bengtsson, Nucl. Phys. A496 (1989) 56. 

[31] G.Andersson, S.E.Larsson, G.Leander, P.Moller, S.G.Nilsson, R.Bengtsson, Nucl. Phys. A268 
(1976) 205. 

[32] K.Neergard, V.V.Pashkevich, and S.Frauendorf, Nucl. Phys. A262(1976) 61. 

[33] S.Aberg, Nucl. Phys. 477 (1988) 18. 

[34] P. Ring and P. Shuck, T/ie Nuclear Many Body Problem (Springer-Verlag, 1980). 

[35] A. Faessler, Fortschr. Phys. 16 (1968) 309. 

[36] I.Hamamoto and H.Sagawa, Nucl.Phys. A327 (1979)99. 

[37] A. Fitzpatrick, S.A. Araddad, R. Chapman, J. Copnell, F. Lind’en, J.C. Lisle, A.G. Smith, J.P. 
Sweeney, D.M. Thompson, W. Urban and S.J. Warburton, J. Simpson, C.W. Beausang, J.F. 
Sharpey-Shafer, S.J. Freeman, S. Leoni, and J. Wrzesinski, Nucl. Phys. A585 (1995)335. 

[38] J.R.B. Oliveira, S. Frauendorf, M.A. Deleplanque, R.M. Diamond, F.S.Stephens, C.W. 
Beausang, J.E.Draper, C.Duyar, E.Rubel, J.A.Becker, E.A. Henry, and N.Roy, Phys. Rev. 
C47(1993)R926. 

[39] T.R. Werner and J. Dudek, Atomic Data and Nuclear Data Tables, 50(1992) 179. 

[40] A. Bohr and B.R. Mottelson, Nuclear Structure, vol. I (Benjamin, 1969). 

[41] R.A. Broglia, T. Dpssing, M. Matsuo, E. Vigezzi, P. Bosetti, A. Bracco, S. Frattini, B. Herskind, 
S. Leoni, and P. Rasmussen, Zeit. Phys. A356 (1996) 259. 

[42] A. Bracco, P. Bosetti, S. Frattini, E. Vigezzi, S. Leoni, T. Dpssing, Herskind, and M. Matsuo, 
Phys. Rev. Lett. 76 (1996) 4484; 

E. Vigezzi, S. Frattini, P. Bosetti, A. Bracco, R.A. Broglia, T. D0ssing,B. Herskind, S. Leoni, M. 
Matsuo, Proc. of the Workshop on Heavy Ion Fusion: Exploring the Variety of Nucl. Properties, 
Padova, 1994, (World Scientific 1995) p.l49. 

[43] M. Matsuo, T. Dpssing, B. Herskind, S. Frauendorf, Nucl. Phys. 564 (1993) 345. 

[44] B.W. Bush, G.F. Bertsch and B.A. Brown, Phys. Rev. C45 (1992) 1709. 

[45] 1. Hamamoto, Nucl. Phys. A232(1974) 445. 

[46] T.Takami, Phys. Rev. Lett. 68 (1992) 3371. 

[47] P.J.Brussard and P.W.M.Glaudemans, Shell-Model Applications in Nuclear Spectroscopy 
(North-Holland, 1977) 


28 



